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ABSTRACT 

We present multi-frequency VLBA 2 observations of PMN J1632-0033, one of 
the few gravitationally lensed quasars suspected of having a central "odd" image. 
The central component has a different spectral index than the two bright quasar 
images. Therefore, either the central component is not a third image, and is 
instead the active nucleus of the lens galaxy, or else it is a third image whose 
spectrum is inverted by free-free absorption in the lens galaxy. In either case, we 
have more constraints on mass models than are usually available for a two-image 
lens, especially when combined with the observed orientations of the radio jets of 
the two bright quasars. If there is no third quasar image, the simplest permitted 
model is a singular isothermal sphere in an external shear field: (3 = 2.05io'io 
(2a), where p(r) oc r _/3 . If the central component is a third image, a hypothesis 
which can be tested with future high-frequency observations, then the density 
distribution is only slightly shallower than isothermal: (3 = 1.91 ±0.02 (2a). We 
also derive limits on the size of a const ant- density core, and the break radius and 
exponent of an inner density cusp. 

Subject headings: gravitational lensing — galaxies: structure — galaxies: nuclei 
- quasars: individual (PMN J1632-0033) 
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1. Introduction 

The central density profile p(r) of galaxies is poorly known, but is important for un- 
derstanding galaxy structure and the nature of dark matter. Theoretical cold-dark-matter 
(CDM) models predict a steep central density cusp, p oc r _/3 with (3 ~ 1-1.5 (see, e.g., 
Navarro, Frenk, & White 1997; Moore 1999). However, measurements of spiral galaxy ro- 
tation curves suggest the dark matter has a constant-density core. This has been argued 
most strongly for low-surface-brightness galaxies (de Blok et al. 2001) but may hold for all 
spiral galaxies (Debattista & Sellwood 2000; Salucci & Burkert 2000), including our own 
Galaxy (Binney & Evans 2001). The conflict may be due to observational effects such as 
beam smearing, at least in part (van den Bosch & Swaters 2001), but the apparent con- 
flict has encouraged some theorists to invent mechanisms that avoid producing cusps (e.g., 
bar-induced evolution, Weinberg & Katz 2002; warm dark matter, Bode, Ostriker, & Turok 
2001; self- interacting dark matter, Spergel & Steinhardt 2000). 

For elliptical galaxies, where rotation curves cannot be measured, the central density 
profile is even less well understood. Observations of early-type galaxies with the Hubble 
Space Telescope (HST) show that the light distribution is cuspy, with cusp parameters that 
may be correlated with global galaxy properties (Lauer et al. 1995; Faber et al. 1997). But 
such studies do not probe mass directly, and are limited to nearby galaxies due to the 
requirements on angular resolution. 

Multiple-image gravitational lenses can be used to measure the central surface density 
of galaxies at significant redshift, by searching for faint central images of the background 
object. In most cases, lenses should produce an odd number of images (Dyer & Roeder 
1980; Burke 1981), one of which is faint and appears near the center of the lens galaxy. 
This "central image" or "odd image" corresponds to the central maximum in the time- 
delay surface. However, of the ~80 known systems in which a galaxy produces multiple 
images of a quasar or radio source, almost all consist of either 2 or 4 images. One system, 
APM 08279+5255, definitely has three quasar images (Lewis et al. 2002b), but the third 
image may be due to a "naked cusp" configuration rather than a central time-delay maximum 
(Lewis et al. 2002a). 

No central image has been identified definitively, even among the radio lenses, where the 
central image could in principle be seen through the light and dust of the lens galaxy. The 
absence of central images has been attributed to the demagnification of the image by the 
high central surface density of the lens galaxy (see, e.g., Narasimha, Subramanian, & Chitre 
1986; Wallington & Narayan 1993; Rusin & Ma 2001; Evans & Hunter 2002; Keeton 2002). 
An "extra" radio component has been detected in several radio lenses, and in three of those 
cases it is still unclear whether the extra component is an additional quasar image or faint 
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radio emission from an active galactic nucleus (AGN). Those three systems are 0957+561 
(Harvanek et al. 1997), MG J1131+0456 (Chen & Hewitt 1993), and the subject of this 
paper: PMN J1632-0033 (Winn et al. 2002). 

Radio components A and B of J1632-0033 are two images of a z = 3.42 quasar, but the 
third, faint component C is of unknown origin (see Fig 1). In an HST image, the lens galaxy- 
appears to be a fairly circular early-type galaxy with an optical effective radius of ~ 0'/2, but 
it is too faint to be characterized in detail. The position of component C is near the optical 
lens galaxy position, as expected for either an AGN or a central odd image. If A, B, and 
C are all images of the same quasar, then they should have similar radio spectra. However, 
Winn et al. (2002) only presented measurements of the flux density of C at one frequency. 
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Fig. 1. — Radio map of J1632-0033 reproduced from Winn et al. (2002), illustrating the con- 
figuration of components A, B, and C. The map is based on the 5.0 GHz merlin observation 
of 2001 April 13. The restoring beam (88 x 44 mas, P.A. 24°) is illustrated in the upper left 
corner. Contours begin at 3a and increase by factors of two, where a = 0.13 mJy beam -1 . 

Here we present new VLB A data that, in combination with the previous data, allow us 
to measure the flux densities of all three components at four widely spaced frequencies. Our 
analysis shows that C has a significantly different radio spectrum than A and B. We argue 
that there are two alternative explanations for this difference. Either C is not a third quasar 
image, and instead is the active nucleus of the lens galaxy; or, C is a third image, but its 
spectral index is inverted by free-free absorption in the lens galaxy. 
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Whichever of these possibilities is true, the additional information from the VLBA data 
provides more constraints on mass models of the lens galaxy than are usually available for 
a two-image lens. If the central component is a third image, then its position and flux are 
valuable constraints. If the central component is the lens galaxy, then our data provide an 
unusually precise position for the lens galaxy, and a strong upper limit on the flux of a third 
image. The VLBA maps also reveal the relative orientations of the radio jets associated with 
A and B. As with all gravitational lensing studies, we still do not have enough information 
to determine the complete density profile of the lens galaxy, but we can use parameterized 
models to explore the implications of the data. 

Historically, two classes of parameterized models have been used to study the central- 
image problem: models with constant-density cores, and models with central cusps. Early 
studies used softened power laws, p oc (r 2 + r 2 ) _/3//2 , and found that central images are 
suppressed if the core radius r c is sufficiently small (see, e.g., Wallington & Narayan 1993). 
Since the demagnification of the central image is determined by the central surface density, 
steeper density profiles (larger (3) allow for larger core radii for the same flux limit on the 
central image. Power laws shallower than isothermal (j3 < 2) produce central images even 
in the limit r c — > 0. 

With the gradual realization that theoretical dark-matter profiles and observed surface- 
brightness profiles of early-type galaxies have cusps rather than cores, attention has shifted 
to the properties of density distributions with central cusps. Rusin & Ma (2001) used the ab- 
sence of odd images in the class sample, presently the largest homogeneous sample of radio 
lenses, to argue that the global mass profiles of the lens galaxies cannot be much shallower 
than isothermal. Muhoz, Kochanek, & Keeton (2001) reached a similar conclusion about 
the particular radio lens CLASS B1933+503. Keeton (2002) found that the distribution of 
stars observed in HST images of nearby early-type galaxies is often sufficiently concentrated 
to suppress the central image without recourse to dark matter. 

We describe our observations of J1632-0033 in § 2, and our interpretation of the results 
in § 3. In § 4, we compute lens models for this system under the two different hypotheses 
for the nature of the central component. We consider density distributions that are scale- 
free power laws, and also determine the required properties of a const ant- density core or an 
inner cusp. Finally, in § 5, we summarize our conclusions and discuss future observations 
that could test more definitively whether PMN J1632-0033 is a three-image quasar. 
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2. Observations 

We observed J1632-0033 with the VLBA on 2002 March 14 and 15, for one eight-hour 
session each day. On the first day, the array included one VLA antenna in place of the VLBA 
antenna at Pie Town, but on the second day, the array consisted of the usual ten antennas. 
During both sessions we alternated between observations at the standard 8.4 GHz (3.6 cm) 
band and the standard 1.7 GHz (18 cm) band. For both bands, the observing bandwidth 
of 32 MHz per polarization was divided into 4 sub-bands. Both senses of polarization were 
recorded with 2-bit sampling. The data were correlated in Socorro, New Mexico, producing 
16 channels of width 500 kHz from each sub-band, with an integration time of one second. 

Calibration was performed with aips 5 using standard procedures. We used the obser- 
vations of J1632-0033 to solve for residual delays, rates, and phases directly (rather than 
using phase-referencing), with a fringe- fitting solution interval of 2 minutes. We used the 
multiple-field clean deconvolution algorithm available in AIPS to deconvolve a field centered 
on A and, simultaneously, a field of the same size centered on the mid-point between B and 
C. For the 8.4 GHz data, the fields were 1024 x 1024 with a scale of 0.2 mas pixel -1 . For 
the 1.7 GHz data, the fields were 512 x 512 with a scale of 1.0 mas pixel -1 . 

At first, the data from each day were analyzed separately. Once we verified that the 
results of the two sessions were consistent, we combined the visibility data from each fre- 
quency to produce final maps. The data were self-calibrated, based on the clean model 
derived from the preliminary maps, with a 15-second solution interval. The final 1.7 GHz 
maps are shown in Fig. 2, and the final 8.4 GHz maps are shown in Fig. 3. In the vicinity of 
components B and C, the noise level in the maps is close to the theoretical limit. However, 
the noise level is three times larger in the vicinity of component A due to large sidelobes. 
Further improvement is difficult because of the near-equatorial location of the object, which 
results in poor coverage of the visibility plane even for an eight-hour Earth-rotation synthesis 
observation. 

The flux densities of A, B, and C at each frequency are reported in Table 1. Also 
included in this table are results based on the previous data of Winn et al. (2002), for those 
cases in which all three components were detected (or could have been detected, given the 
angular resolution). These previous data were re-analyzed using natural weighting rather 
than uniform weighting, to maximize the sensitivity to the faint component C and to match 
the present analysis. The results were consistent with previous results in all cases. Of 



5 The Astronomical Image Processing System (aips) is developed and distributed by the National Radio 
Astronomy Observatory. 
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Fig. 2.— VLBA maps of J1632-0033 at 1.7 GHz. The restoring beam (14.3 x 6.1 mas, 
P.A. 0?4) is illustrated in the lower left corner of each panel. Contours begin at 3a and 
increase by factors of two. The dashed contour is — 3a. In the left panel (components B and 
C), a — 0.045 mJy beam -1 . In the right panel (component A), a — 0.135 mJy beam -1 . The 
dotted lines are the la limits we adopted for the jet position angles. 
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Fig. 3.— VLBA maps of J1632-0033 at 8.4 GHz. The restoring beam (2.8 x 1.2 mas, 
P.A. — 0?5) is illustrated in the lower left corner of each panel. Contours begin at 3<r and 
increase by factors of two. The dashed contour is —3a. In the left and center panels 
(components B and C), a = 0.040 mJy beam -1 . In the right panel (component A), a = 
0.120 mJy beam -1 . The dotted lines are the la limits we adopted for the jet position angles. 
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particular interest was the improved 15 GHz map, in which component C was detected with 
flux density 0.84 ± 0.19 mJy. This is consistent with the upper limit of 1.5 mJy given by 
Winn et al. (2002), but the positive detection of C in the lower-noise map provides another 
valuable point of comparison. 

For the VLA and merlin 6 data, the flux densities were measured by fitting a model 
consisting of 3 point sources to the visibility data, using difmap (Shepherd, Pearson, & 
Taylor 1994). The VLBA data sets were too large for model-fitting in visibility space. For 
these data, the flux densities were measured in image space with aips, using an aperture 
surrounding each component. The quoted uncertainties are the larger of the rms noise and 
the spread in values obtained using different choices for the aperture. 

With these measurements of the flux densities of all three components at four widely 
spaced frequencies, we have the opportunity to test whether C has the same radio continuum 
spectrum as A and B. Assuming that the flux density of each component obeys Si(u) oc uf 
(where a, is the spectral index of component i), a linear least-squares fit to the data in Table 1 
gives a A = -0.27 ± 0.07, a B = -0.25 ± 0.07, and a c = +0.29 ± 0.18. The uncertainties in 
the spectral indices include a contribution from the uncertainties in the absolute flux density 
scale, which is at least 5% at all frequencies. 

The comparison between cka and ac can be made more precise by examining the 
flux density ratios, because the ratios are independent of the uncertainties in the abso- 
lute flux density scale. The difference in spectral indices, ac — «a, is the logarithmic 
slope of Sc(v)/Sa(v)- The flux density ratios Sb(v)/Sa(v) and Sc^V-Sa^) are plotted 
in Fig. 4. The solid line represents the best-fit power law to Sc(v)/Sa(v), which has a slope 



6 The Multi-Element Radio Linked Interferometry Network (merlin) is a U.K. national facility operated 
by the University of Manchester at Jodrcll Bank Observatory on behalf of PPARC. 



Table 1. Flux density measurements of J1632-0033 



Date 


Observatory 


Frequency 




Flux density (mJy) 






(GHz) 


A 


B 


C 


2002 Mar 14-15 


VLBA 


1.67 


242 ±3 


16.9 ±0.1 


0.45 ±0.05 


2000 Apr 29 


VLBA 


4.99 


167 ±2 


11.9 ±0.1 


0.60 ±0.11 


2001 Apr 05 


MERLIN 


4.99 


207 ±2 


14.2 ± 0.3 


0.90 ±0.11 


2002 Mar 14-15 


VLBA 


8.42 


167 ±2 


14.0 ±0.3 


0.73 ±0.04 


2000 Nov 11 


VLA 


14.9 


144 ± 2 


9.4 ±0.2 


0.84 ±0.19 


2000 Nov 11 


VLA 


22.5 


126 ± 2 


8.4 ±0.6 


< 3.2 


2000 Nov 11 


VLA 


42.6 


83 ± 1 


5.4 ±0.4 


< 1.8 
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ac — «a = +0.51 ± 0.17, implying the spectral indices of C and A are discrepant at the 
3a level. The dashed line is a fit to the data under the hypothesis that C is a third quasar 
image whose spectrum had been modified by free-free absorption (see § 3). Also plotted are 
upper limits on Sq(v)/ Sa(v) from VLA data at the two highest frequencies. 

The position angles of the jets in A and B were measured by fitting elliptical components 
to the VLBA maps with aips, using a procedure (JMFIT) that corrects for the ellipticity 
of the beam. For modeling purposes, we adopted the average of the 1.7 GHz and 8.4 GHz 
results, giving a — 8° and 0b — —69°, measured east of north. We estimate the uncertainty 
in each position angle to be ±15°, an error range that is approximately twice the difference 
between the 1.7 GHz and 8.4 GHz results. These ranges of position angle are illustrated by 
dotted lines in Figs. 2 and 3, although it is difficult to account for the elliptical beam in a 
visual inspection. 

3. The nature of component C 

The flux density ratio Sb^/Sa^) is nearly independent of observing frequency, as 
should be the case for gravitationally lensed images of a single compact source. By contrast, 
the flux density ratio Sciy)/ Sa{v) increases significantly with frequency. This immediately 
suggests that C is not a third image of the same background source. Given that its position is 
consistent with the optical position of the lens galaxy, it would be natural to conclude that C 
represents faint radio emission from the lens galaxy, as discussed further in § 3.1. However, 
there are several mechanisms that might cause a third quasar image to appear to have a 
different spectral index than the two brighter images. These are discussed in § 3.2— § 3.5. 
The only mechanism we cannot rule out is free-free absorption (§ 3.5). 

3.1. Lens galaxy emission 

How does the flux density of C compare with the flux density expected from a typical 
galaxy at the lens redshift? Given the radio flux density and spectral index of component 
C that were presented in § 2, and the lens redshift estimate 7 of z ~ 1, the implied 1.4 GHz 
luminosity of the lens galaxy is 1.7 x 10 24 /i~ 2 W Hz" 1 . 



7 The lens redshift zi = 1.0 ± 0.1 was estimated by Winn et al. (2002) by requiring the photometric 
properties of the lens galaxy to conform to the fundamental plane of elliptical galaxies, a method invented 
by Kochanek et al. (2000). 
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Fig. 4. — Flux density ratios Sb(v)/Sa{v) and Sc(v)/Sa(is), as a function of frequency. To 
allow both quantities to be plotted on the same axes, the values of Sc(v)/Sa(is) have been 
multiplied by 10. Upper limits on Sci^/S^u) are plotted at 22.5 GHz and 43 GHz. The 
straight line is the best-fitting power law to the Sq{u) / Sx{u) data points at 15 GHz and 
below. The dotted line is the best-fitting curve assuming the variation is entirely due to 
free-free absorption (see Eq. 1 and § 3.5). 
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This is significantly more powerful than the radio emission observed from ordinary galax- 
ies, which is thought to arise from star formation and is limited to £ 10 23 W Hz -1 . However, 
it is within the observed range of luminosities from galaxies harboring active galactic nuclei 
(AGN). For example, Sadler et al. (2002) identified a sample of 420 AGN in the 2dF Galaxy 
Redshift Survey, of which 24% had a 1.4 GHz radio luminosity at least as large as the implied 
luminosity of component C (after correcting for their different choice of cosmology). 

We conclude that the identification of component C as the active galactic nucleus of 
the lens galaxy is a reasonable one. This would make J1632-0033 unusual among the known 
lenses, but not unique; the lens galaxy of B2045+265 has a bright central component that 
is almost certainly an active nucleus (Fassnacht et al. 1999). 

3.2. Variability 

If the source is variable, then the instantaneous flux density ratios between the lensed 
images will fluctuate, due to the differential time delays between the images. If, in addition, 
the degree of variability depends on frequency, then the flux density ratios will vary with 
frequency. However, the differential time delay between components B and C is expected to 
be very small (<0.5 day). This is much smaller than typical time scales for quasar variability 
(~months). Indeed, our measurements, which were spaced apart by one day, showed no sign 
of variability. We therefore discount the possibility that variability is the reason for the 
discrepancy of spectral indices. 

3.3. Radio substructure 

If the background radio source is extended and has sub-components with different spec- 
tral indices, then these sub-components may each be magnified by a different factor (see, 
e.g., Patnaik & Porcas 1999). In that case, the sum of flux densities of the sub-components 
in a lensed image would not necessarily have the same spectral index as the corresponding 
sum in a different image. 

However, in the case of J1632-0033, the source is very compact. From Fig. 3 we can 
limit the separation of sub-components in image A to <5 mas, corresponding to a change in 
magnification of only A/i ^ 0.01 in an isothermal model. Furthermore, the agreement of «a 
and « B argues against the presence of a significant effect from frequency-dependent radio 
substructure. 
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3.4. Resolved-out radio structures 

If a radio source has structure on large angular scales, with Fourier components corre- 
sponding to smaller spatial frequencies than those measured by an interferometer, then that 
portion of the radio source is invisible to the interferometer. The invisible radio structure 
is said to be "resolved out." We must consider the possibility that the spectral indices of A 
and C are different only because the components are being resolved out by different amounts 
at each frequency. 

We reject this hypothesis based on three facts. First, there cannot be much resolved- 
out structure in components A and B. During a monitoring program, we measured their 
flux densities at 8.4 GHz with the VLA on 2002 March 14, simultaneous with the VLBA 
measurement described in § 2. The VLBA flux densities were 83 ± 7% of the VLA flux 
densities, despite being measured on angular scales ~200 times smaller. Second, if C is a 
third quasar image, it is demagnified by a factor of ~10 2 and would have an even smaller 
fraction of resolved-out structure than A or B. Third, even if there is resolved-out structure 
in A, it would be expected to have a steep radio spectrum, because the extended jets of 
quasars generally have steeper spectra than the compact cores. This would only worsen the 
discrepancy between the spectral indices of A and C, by causing the true spectral index of 
A to be steeper than observed. 

3.5. Propagation effects 

Propagation effects that are chromatic, and that affect each lensed image by a different 
amount, will cause the spectra of the images to differ. For example, reddening by dust 
causes the optical colors of lensed images to differ (see, e.g., Motta et al. 2002; Falco et al. 
1999). At radio wavelengths, the propagation effects are due to interstellar plasma. The 
polarizations of lensed images may differ due to Faraday rotation (see, e.g., Patnaik et al. 
1993). Or, in a phenomenon known as "scatter-broadening," the angular size of a lensed 
image may be increased due to scattering by electron-density fluctuations along the line of 
sight (see, e.g. Jones et al. 1996). However, these two effects do not alter the spectral indices 
of lensed images. A possible exception might arise from scatter-broadening, if it causes a 
frequency-dependent amount of flux to be resolved out (see § 3.4 above), but our data do 
not show the characteristic v~ 2 dependence of angular size due to scatter-broadening. 

One propagation effect that would alter the spectral indices of lensed images is free- 
free absorption by electrons in the lens galaxy. Free-free absorption decreases strongly with 
frequency, which would flatten or invert the spectral index of the central component. We 
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must therefore consider the possibility that component C is a third quasar image whose 
spectral index has been inverted by free-free absorption. 

Assuming that C is affected by free- free absorption, and A is not, we would expect the 
flux density ratio to obey Sc(v)/Sa(v) = /iqa exp(— t u ), where /xqa is the lensing magnifi- 
cation ratio and r„ is the optical depth to free-free absorption. Here v — (1 + z)u Q \, s is the 
frequency at the lens redshift. A useful approximation for r„ has been given by Altenhoff et 
al. (1960) and Mezger & Henderson (1967): 

/T e V L35 / v \-2.i/ [N?ds\ (v\~ 21 
r„ = 0.08235 x(^j (^) {^- 6 ) = {-) , (1) 

where T e and iV e are the electron temperature and number density, and v c is defined im- 
plicitly. This formula predicts a sharp cut-off at low frequencies, and the asymptotic limit 
Sc{v)lS K {v) = fi CA as v — > oo. 

We solved for the values of the two parameters u c and /icA that provide the best fit to 
the data. The results were v c = 3.2 GHz (z/ Ci0bs = 1.6 GHz) and /j,ca = 0.0045. The dotted 
line in Fig 4 shows the best-fitting curve, which is statistically consistent with the data 
(x 2 = 11 with 3 degrees of freedom). There are no published cases of free- free absorption 
by a lens galaxy for comparison. However, the best-fitting value v c = 3.2 GHz is comparable 
to the values of v c = 1-10 GHz that have been inferred for the central regions of some active 
galaxies, using multi-frequency VLBI measurements of radio jets (for recent examples, see 
Marr, Taylor, & Crawford 2001, Jones et al. 2001, Kameno et al. 2000, Taylor 1996, and 
Levinson, Laor, & Vermuelen 1995). 

Because the data are quantitatively consistent with both the third-image hypothesis and 
the free-free absorption hypothesis, we have only our prior prejudice on the relative likelihood 
of these two possibilities to guide us. One might object that the third-image hypothesis 
is contrived, because neither a central image nor free-free absorption has previously been 
observed in a radio lens. However, we do not consider this "double novelty" as a serious 
objection because a central image might be expected to suffer more absorption than other 
images, due to the potentially large electron column density near the center of the lens 
galaxy. For this reason, we take a neutral position, and consider the consequences of both 
hypotheses in the section that follows. 

4. Updated lens models 

In this section we use information from the VLBA maps to update the lens models 
for J1632-0033 that were presented by Winn et al. (2002). The goal of our analysis is to 
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investigate the improvement that results when two extra pieces of information are available 
besides the usual image positions and fluxes: the relative orientations of the radio jets, and 
a strong constraint (either a measurement or an upper limit) on the flux of a central image. 
In § 4.1, we assume that component C is due to the lens galaxy. In § 4.2, we explore the 
consequences if component C is a third quasar image. 



4.1. Case 1: Component C is the lens galaxy 

In this case, we use the model constraints given in Table 2. The position of C is taken 
to be the center of the lens galaxy, which is known more precisely than the lens galaxy 
position in most other systems. The upper limit on the flux density of a third quasar im- 
age, expressed as the magnification ratio fi^/fix, represents the 5a limit in the 1.7 GHz 
VLBA map. The magnification ratio /iba = He/ Ha is the average of the 15 measurements 
of Sb(v)/ Sa(v) presented in this paper and by Winn et al. (2002), and the uncertainty is 
the variance in those measurements. In theory, our error estimates for the positions and 
magnification ratio should be enlarged (to ~l-2 mas and 10-20%, respectively) to account 
for possible mass substructure in the lens galaxy (see, e.g., Metcalf & Madau 2001; Dalai 
& Kochanek 2002). However, in all the cases described below, the error statistic is domi- 
nated by either the third-image flux or the jet orientations, causing the precise values of the 
uncertainties on the positions and magnification ratio to be unimportant. 



We define the goodness-of-fit parameter as 



* 2 = £ 



i=A,B 



(AR,)^ ' (A4>i) 



+ 



, (/fgA ~ A^Ba) 2 / 9 x 
+ (A/iBA) 2 ' {) 



in which primed quantities are those determined by the model, and unprimed quantities are 
those observed. In addition, we reject models that predict a third image brighter than the 
upper limit. This is computationally faster than the more orthodox procedure of adding 
another error term to Eq. 2 for the third-image flux, but we have verified that the results 
are very nearly equivalent. This is because the flux of the third image depends sensitively 
on the parameters in our models, causing x 2 to increase very rapidly across the boundary in 
parameter space where the third-image constraint is violated. 

Although in principle there could be a complicated angular structure in the lens model, 
due to the ellipticity of the lens galaxy and due to external perturbations, we do not have 
enough constraints to explore such models. Instead, we use mass models with circular 
symmetry, to which we add an external shear field (of strength 7 and position angle # 7 ) 
to simulate the combined effect of galaxy ellipticity and tidal fields from neighboring mass 
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concentrations. 

We begin with global power-law models, for which p(r) oc r _/3 . The model has 7 free 
parameters: (3, 7, 9 y , the mass normalization, the source location (x s ,y s ), and the position 
angle of the source jet <f> B . We applied the 8 constraints of Table 2 using a ^-minimization 
code. Because the constraint on the third-image flux is one-sided, and because of the way in 
which we have implemented that constraint, the allowed models fit the constraints perfectly 
(x 2 = 0). The optimal value is f3 = 2.05 and the 2a bounds are 1.95 < j3 < 2.28. This 
range brackets the isothermal value (3 = 2. The upper bound is enforced by the observed jet 
orientations. The lower bound is enforced by the upper limit on the flux of the third image. 

The lower bound is robust, even if our upper limit on // 3 / /j, a is too strong due to the 
possible extinction of a third image by free-free absorption, or the lowering of its peak flux 
density by scatter-broadening. For example, if we weaken the limit on the third-image flux 
density by a factor of 10 (//3///A < 0.01), the lower bound changes only slightly: 1.91 < (3 < 
2.28. 

Interestingly, there is a different (weaker) lower bound on f3 that results from the relative 
positions of A and B with respect to the lens center. For f3 < 1.83, the radial critical curve is 
sufficiently large to encompass B. In such models, the parity of B is reversed, and a brighter 
third image is predicted in a location outside the radial critical curve; image B becomes the 
central odd image in a three-image system. Such models can be ruled out. 

Next, we investigate the possibility of a const ant- density core. Given the preceding 
results, we assume that the density profile at large radii is isothermal and derive an upper 
bound on the size of a hypothetical const ant- density core. We adopt a two-dimensional 
surface density of the form 

k(R) = . b (3) 

2y/WTRj 

where R is the two-dimensional coordinate. The constraint on the third-image flux requires 
R c < 1.6 mas, corresponding to 0.8% of the optical radius of the lens galaxy, or 9k" 1 parsecs 
at z = l. 2 

Finally, we explore models with central cusps. We use the cusped models introduced by 
Muhoz, Kochanek, & Keeton (2001), in which the mass density is 

Pij) = (r/r b )P[l + (r/r b ) 2 ]b-W (4) 



2 Here and elsewhere, we compute angular-diameter distances assuming a flat fl m — 0.3 cosmology, giving 
5.61/i -1 kpc per arc second. 
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The density falls off with a logarithmic slope of rj for r 3> r&, and /3 for r <C r^, where r?, is 
the "break radius." This model is convenient because the lensing deflection can be computed 
analytically. 

First, we assume that the mass distribution at large radii is isothermal (rj — 2). We use 
the constraints in Table 2 to find the allowed region in the two-dimensional space spanned 
by (3 and r&. As before, the allowed models fit the constraints perfectly (x 2 = 0). The 
results are shown in the left panel of Fig. 5. The la allowed region is shaded. For large 
Tf,, the overall mass distribution must be nearly isothermal, as expected from the preceding 
results. Shallower inner cusps are allowed only when ^ 0'.'02, corresponding to 10% of the 
optical radius, or 100/i -1 parsecs at z = 1. As r& approaches zero, the constraint on f3 widens 
rapidly. At r b m 1.6 mas the range of allowed exponents encompasses (3 = 0, which is the 
limit of a constant-density core derived in the previous section. 

The allowed region of cusp parameters is qualitatively different if we allow the mass 
profile of the galaxy at large radii to be steeper than isothermal. The upper and lower right 
panels of Fig. 5 show the results for the cases rj = 3 and rj = 4. Note that both NFW 
(Navarro, Frenk, & White 1997) and Hernquist (1990) profiles, with (3 = 1 central cusps and 
asymptotic slopes of rj = 3 and rj = 4 respectively, are not permitted as global models of the 
mass distribution, because their inner cusps are too shallow. Unlike the case of rj = 2, the 
range permitted for the central cusp exponent j3 does not widen as r& shrinks. The central 
image and jet orientation constraints instead pinch off the permitted parameter space to set 
a minimum break radius: r& > 0'/12 (rj, > 0'.'2) for rj = 3 (77 = 4). Furthermore, for both 
rj = 3 and 77 = 4, the allowed values of f3 bracket the isothermal value (3 — 2, regardless 
of the break radius. The non-isothermal profile at large radii must be compensated by a 
nearly isothermal profile within a fairly large break radius. Again, we see that if the mass 
profile obeys a power law over the range of radii between images A and B, then the favored 
distribution is close to isothermal. 



4.2. Case 2: Component C is a third image 

For the three-image scenario, we use the model constraints given in Table 3. In this 
case, the lens galaxy position is taken from the optical measurement of Winn et al. (2002), 
with error bars enlarged to 23 mas (0.5 pixel) to allow for systematic errors. 3 The value 



3 The smaller error bars of 9 mas quoted by Winn et al. (2002) assume a particular form for the point- 
spread function and the galaxy profile. Ordinarily, this would not introduce much additional error, but in 
this case component B and the lens galaxy are both faint and merged in the HST image. 
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Inner power-law slope /S Inner power-law slope /J 



Fig. 5. — Cuspy lens models, assuming that component C is the lens galaxy. In all cases, 
the la allowed region is shaded. The left boundary of the shaded region represents the 
third-image constraint. The right boundary of the shaded region is the la boundary due 
the jet orientations (the 2a boundary is also shown). The dashed curve is where image B 
flips parity. The horizontal dotted lines mark the observed radii of images A and B, and the 
/-band effective radius of the lens galaxy measured by Winn et al. (2002). 



adopted for Sc(v)/Sa(is) is based on the fit to the free-free absorption model described in 
§ 3.5. We define the goodness-of-fit parameter as 



x 2 = E 



i=B,C,Gal 



~ (r:-r.) 2 

(ARj) 2 



i=A,B L 



(0; 




+ 



(a^ba ~ /^ba) 2 (/^ca ~ /^ca) 2 
(A/i BA ) 2 ~ (A/i CA ) 2 



(5) 



in which primed quantities are those determined by the model, and unprimed quantities are 
those observed. 

For a scale-free power-law profile in an external shear field, we have three degrees of 
freedom. The model does a good job of reproducing the configuration, with \ 2 — 2.2. The 
logarithmic slope of the density profile in the optimal model is f3 — 1.91 ± 0.02 (2a). If we 
take the radial density profile to be isothermal with a constant-density core, as in Eq. 3, we 
can also achieve a reasonable fit (x 2 = 3.0), with the optimal value r c = 5.5 milliarcseconds, 
or 31/i -1 parsecs, and a 2a confidence region of 3.6-7.8 mas. 

When we apply the cuspy lens model of Eq. 4, we find that the allowed region of cusp 
parameters is much more tightly constrained than they are in the two-image lens models. 
Figure 6 shows the results for models in which the radial density distribution at large radii 
is isothermal (left panel), or obeys a power law with logarithmic slope rj = 3 (upper right), 
or r] = 4 (lower right). Both boundaries of the allowed region (shaded) are enforced by the 
constraint on the third-image flux. In particular, for (3 > 2, no third image is produced. 
For large r b , the cusp must conform to the value (3 = 1.91 derived above. As r b becomes 
smaller, the inner mass distribution must become shallower to compensate for the steeper 
mass distribution at large radii. 



For the purpose of finding central odd images, the best systems are radio-loud two-image 
quasars with a large magnification ratio between the two bright images. Radio loudness is 
important because the radio components can be observed through the lens galaxy, whereas 
optical components are likely to be hidden by the lens galaxy. Free-free absorption is poten- 
tially significant but can be overcome by observing at high frequencies. Two-image systems 
are better than four- image systems, as emphasized recently by Keeton (2002). This is be- 
cause the projected source position in four-image systems is usually closer to the lens center, 
where the central image magnification is smallest. In addition, four-image systems are sub- 
ject to a larger magnification bias than two-image systems (Wallington & Narayan 1993; 
King & Browne 1996; Keeton, Kochanek, & Seljak 1997; Rusin & Tegmark 2001; Finch et 



5. Discussion and summary of conclusions 
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Fig. 6. — Cuspy lens models, assuming that component C is a third quasar image. In all 
cases, the la allowed region is shaded. The la and 2a boundaries are marked by solid lines. 
The horizontal dotted lines mark the observed radii of images A and B, and the J-band 
effective radius of the lens galaxy measured by Winn et al. (2002). 
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al. 2002), leading to a smaller intrinsic source flux and a fainter odd image. Among two-image 
systems, the asymmetric doubles are favorable because the source position is close to the 
radial caustic, producing the brightest possible central image for a given mass distribution. 

The subject of this paper, PMN J1632-0033, is exactly this type of system: an asym- 
metric two-image radio lens. For this reason, it seemed reasonable that the third component 
found by Winn et al. (2002) could be an example of the elusive and long-sought central odd 
images. However, we have shown that the third component has a spectral index that differs 
by 3<t from those of the two bright quasar images. 

We have argued that the two most plausible interpretations are that the central com- 
ponent is emission from the lens galaxy, or else that it is a third image with the extra 
complication that its spectral index has been modified by free-free absorption. The data are 
quantitatively consistent with either scenario. Lens galaxies are rarely radio-loud, but the 
flux density of the component is consistent with typical radio powers of active galactic nuclei. 
No central image has ever been securely identified for a radio lens, nor has free-free absorp- 
tion by a lens galaxy ever been reported (with one possible exception, class B0128+437; 
I. Browne 2002, private communication), but one might expect the two phenomena to be 
related due to the hot and dense conditions near galaxy centers. 

Given the new VLBI data, we have investigated models for the mass distribution of the 
lens galaxy under various assumptions. Whether or not component C is a third image, we 
showed that if the mass distribution obeys a power law p oc r _/3 , then it must be nearly 
isothermal. If C is the lens galaxy then 1.95 < (3 < 2.28. If C is a third image then the 
profile is tightly constrained to be only slightly shallower: 1.89 < /3 < 1.93. We also derived 
limits on a constant-density core radius, or on the break radius and inner exponent of a 
central cusp. 

These results can be added to the growing collection of evidence that early-type galaxies 
have nearly isothermal mass profiles. In the few other gravitational lens systems where it is 
possible to measure the radial density profile, the results favor isothermal profiles (Kochanek 
1995; Cohn et al. 2001; Muhoz, Kochanek, & Keeton 2001; Rusin et al. 2002; Koopmans 
& Treu 2002), although there are some possible counter-examples (Treu & Koopmans 2002; 
Chen, Kochanek, & Hewitt 1995). Moreover, the findings are consistent with studies of 
early-type galaxies based on stellar dynamics (see, e.g., Gerhard et al. 2001; Rix et al. 1997) 
and X-ray halos (see, e.g., Fabbiano et al. 1989). 

To further elucidate the nature of the central component and the central density profile 
of the lens galaxy, it would help to have a better optical image of the lens galaxy. The 
current HST image provides only a weak detection of the lens galaxy. A better image could 
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be used to test more precisely whether radio component C is coincident with the center of 
the lens galaxy, and to study its surface brightness profile in detail. For example, with a 
better measurement of the galaxy light profile, one might be able to explore the balance 
between the luminous and dark matter in this system. 

Even more helpful would be sensitive observations at high radio frequencies ( £ 30 GHz). 
As shown in Fig. 4, the power-law and the free-free absorption fits diverge significantly at 
high frequencies, where free-free absorption becomes negligible. For free-free absorption, 
the ratio Sc(v) / Sa(v) must converge to a constant value. For lens galaxy emission, the 
spectrum need not follow a power law at the high frequencies, but it would require an 
unlikely coincidence for the intrinsic radio spectrum of the lens galaxy to turn over at high 
frequencies and become identical to the radio spectra of A and B. Observations at low radio 
frequencies ( £ 1 GHz) could also distinguish the hypotheses, in principle, but it would be 
harder to achieve the necessary angular resolution and signal-to-noise ratio. High frequency 
observations present their own challenges, including high noise levels and susceptibility to 
poor atmospheric conditions, but the reward would be a definitive discovery of one of the 
long-sought but never-seen central images of a gravitational lens. 

We thank the anonymous referee for constructive criticism. We also thank I. Browne and 
L. Koopmans for interesting discussions about radio propagation effects. J.N.W. is supported 
by an Astronomy & Astrophysics Postdoctoral Fellowship, under NSF grant AST-0 104347. 
C.S.K. is supported by NASA ATP grant NAG5-9265. 
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Table 2. Constraints on lens models in which C is the lens galaxy 



Parameter 


Value 


R.A.(B) - R.A.(A) 


+1241.90 ±0.07 mas 


Decl.(B) - Dccl.(A) 


-784.00 ±0.16 mas 


R.A.(Gal) - R.A.(A) 


+1160.61 ±0.07 mas 


Decl.(Gal) - Decl.(A) 


-726.10 + 0.16 mas 


Hba = Mb/ma 


0.076 + 0.008 


M3/MA 


< 0.00093 


4>a = P A. of A-jct 


8 + 15° 


4> B = P.A. of B-jet 


-69 ± 15° 



Table 3. Constraints on lens models in which C is a third image 



Parameter 


Value 


R.A.(B) - R.A.(A) 


+1241.90 + 0.07 mas 


Decl.(B) - Decl.(A) 


-784.00 ± 0.16 mas 


R.A.(C) - R.A.(A) 


+1160.61 ± 0.07 mas 


Decl.(C) - Decl.(A) 


-726.10 ±0.16 mas 


R.A.(Gal) - R.A.(A) 


+1162 ± 23 mas 


Decl.(Gal) - Decl.(A) 


-738 ± 23 mas 


Hba = Mb/ma 


0.076 ±0.008 


MCA = Mc/^a 


0.0045 ± 0.0009 


4>a = P.A. of A-jct 


8± 15° 


4> B = P.A. of B-jet 


-69 ± 15° 



